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ABSTRACT 

We present the UDF05 project, a HST Large Program of deep ACS (F606W, F775W, F850LP) and 
NICMOS (FllOW, F160W) imaging of three fields, two of which coincide with the NICPl-4 NICMOS 
paraUel observations of the Hubble Ultra Deep Field (HUDF). In this first paper we use the ACS 
data for the NICP12 field, as well as the original HUDF ACS data, to measure the UV Luminosity 
Function (LF) of z ~ 5 Lyman Break Galaxies (LBGs) down to very faint levels. Specifically, based 
onaV — i, i — z selection criterion, we identify a sample of 101 and 133 candidate z 5 galaxies down 
to Z850 = 28.5 and 29.25 magnitudes in the NICP12 and in the HUDF fields, respectively. Using an 
extensive set of Monte Carlo simulations we derive corrections for observational biases and selection 
effects, and construct the rest-frame 1400 A LBG LF over the range M1400 = [—21.4,-17.1], i.e. 
down to ^ 0.04 L^, a t z ~ 5, and complement it with data from the Subaru Deep Field (SDF) from 
lYoshida et"aH ()2006l ) to extend it to the brighter end (M1400 > -22.2). We show that: (i) Different 
assumptions regarding the SED distribution of the LBG population, dust properties and intergalactic 
absorption result in a 25% variation in the number density of LBGs at 2; ~ 5; (ii) Under consistent 
assumptions for dust properties and intergalactic absorption, the HUDF is about 30% under-dense 
in z ~ 5 LBGs relative to the NICP12 field, a variation which is well explained by cosmic variance; 
(iii) The faint-end slope of the LF is independent of the specific assumptions for the input physical 
parameters, and has a value of a ~ —1.6, similar to the faint-end slope of the LF that has been 
measured for LBGs at z ~ 3 and z Q. Our study therefore supports no variation in the faint-end of 
the LBG LF over the whole redshift range z ~ 3 to z '--^ 6. Based on a comparison with semi-analytical 
models, we speculate that the z ~ 5 LBGs might have a top-heavy IMF. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
luminosity function — dark matter 



I. INTRODUCTION 

Deep imaging with the Advanced Camera for Surveys 
(ACS) and the NICMOS camera onboard the Hubble 
Space Telescope (HST) has been instrumental in push- 
ing the study of galaxy populations out to the reion- 
ization frontier of z ~ 6 and beyond. The develop- 
ment of the Lyman Break technique has been a mile- 
stone for the study of high redshift galaxies, as it al- 
lows to identify, from broad-band photometry alone, 
large numbers of star-forming, but otherwise normal, 
galaxies at early epochs: star forming galaxies at a 
given high redshift are identified from the drop in fiux 
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by, typically, 1-2 magnitudes, blueward of the detec- 
tion passband, caused by the Lyman continuum break 
and the Lyman series blanketing induced by intergalac- 
tic neutral hydrogen clouds. While th e conceptual 
idea was worked out already 30 years ago (|Meier. D. L.I 
|1976[) . it was only in the 1990's that progress in in- 
strumentation made it possible t o identify high-z galax- 
ies using this techiiique ( e.g. ICuhathakurta et al.lll990t 
ISteidel k, Hamiltoiilll992[ ). Several spectroscopic follow- 
on studies of LBG candidates have proven the high ef- 
ficiency and reliability of the Lyman Break tec hnique 
in identifying galaxies at early epochs (e.g, . Steidc l et al.l 
199^; iMalhotra et ani2005l : iVanzella et aLlfeOOQ) . 

Despite the availability of large samples of high red- 
shift galaxies with low contamination from interlopers, 
the derivation of a Luminosity Function (LF), and par- 
ticularly of its faint-end slope, remains difficult on ac- 
count of the small volume probed by currently available 
ultra de ep surveys such as th e Hubble Ultra Deep Field 
(HUDF; iBeckwith et"aIll2006D . and because of the many 
corrections which are often not well constrained. This 
is illustrated by the large scatter in the published LF 
parameter s that have been derived from the same data 
sets fe.g. iBunke r et al.l 120 041: lYan fc Windhorsti 12004 
iBouwens et al.ll2 006: BcckwitheLalJl2006D. Mapping the 
evolution of the faint-end slope of the LF at high redshifts 
remains a major goal of observational cosmology, as the 
faint-end slope is expected to be dramatically affected 
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Fig. 1. — The planned UDF05 fields relative to the origi- 
nal HUDF ACS and NICMOS-parallel pointings. The larger and 
smaller squares correspond to the ACS and NICMOS pointings, 
respectively. The Veoe -dropout study that we present in this pa- 
per is based on the UDF05 ACS observations covering the NICP12 
field, and the HUDF. Due to the two-gyro operation of the HST, 
the UDF05 adds a new deep-imaging NICMOS field to the original 
HUDF pointings (NICP5) . 

by reionization (e.g. IWvithe fc Loebll2006[ ). and a solid 
detection of a change in slope would be strong evidence 
for the transition between a partially and a fully reion- 
izcd universe. As also highlighted by e.g., Bcckwith et_alj 
([2OO6), one important step to improve the sensitivity to 
changes in the faint end slope is to derive LFs at different 
redshifts in the most homogeneous way, applying similar 
techniques to similar data sets. 

This paper presents the UDF05 project, a 204-orbit 
HST Large Program of ultradeep ACS and NICMOS ob- 
servations of multiple fields to study the evolution of the 
faint-end of the LF throughout the z A — 8 redshift 
regime. The UDF05 was originally constructed to ob- 
serve with ACS/WFC the two NICMOS parallel fields 
that were acquired while the HUDF was imaged with 
the ACS. These HUDF-NICMOS parallel fields (here- 
after NICP12 and NICP34, see Fig. [U, have a total ex- 
posure of ~ 2 X 10^ s both in the FllOW (broad J) and 
F160W (H) filters, providing the deepest ncar-IR data 
available. Only relatively shallow ACS data were avail- 
able at their locations, and our additional UDF05 data 
for these fields was aimed at reaching a depth in the visi- 
ble, comparable to that of the NICMOS observations (see 
Table 1). Furthermore, the planned UDF05 pointings 
should have had an orientation such that, while acquir- 
ing the ACS data for the NICP12 and NICP34 fields, 
parallel NICMOS observations should have been taken 
for the original HUDF, so as to obtain a total of three 
fields with ultradeep imaging both in the visible and in 
the near-IR. Unfortunately, the transition of HST to the 
two-gyro mode and the first failure of ACS in June 2006 
severely limited the number of available orbits at the re- 
quired orientation, so that the NICMOS observations of 
the original ACS HUDF field ended up with only one 
third of the exposure time originally planned. 

In our first utilization of the UDF05 dataset we study 
the faint-end of the LF of z ~ 5 Vf^oe-dropout galaxies; 
specifically, we compare the z ~ 5 LF obtained from the 
ACS NICP12 field with that derived from the original 
ACS HUDF, in order to assess the effects of cosmic vari- 
ance on the z ^ 5 LF parameters. In Oesch et al. (2007, 
in preparation) and Stiavelli et al. (2007, in preparation) 



we will extend our exploration of the realistic uncertain- 
ties in the faint-end slope of the LF of Lyman Break 
Galaxies to the z ~ 6 and z > 7 regimes, respectively, 
probing well into the expected reionization epoch and 
thus constraining its impact on galaxy formation. 

This paper is structured as follows. After the descrip- 
tion of the data, corrections for noise correlation and 
source detection (^J2]), we describe our selection criteria 
(iSJ, and the simulations that were performed to esti- 
mate the completeness of our source catalogs O. The 
observed surface densities are derived by including the ef- 
fects of photometric errors (^ and the LF is computed 
under different assumptions about the underlying SED 
distribution of the LEG population as well as different 
prescriptions for intergalactic hydrogen absorption O. 
In section f}7] we present a theoretical estimate for the 
effects of cosmic variance on our results. In section ^ 
we discuss the evolution of the faint-end of the LF of 
star forming galaxies in the z ~ 3 — 6 redshift regime, 
and compare our observed z ~ 5 LF with theoretical 
predictions based on a standard^" A Cold Dark Matter 
(A-CDM) universe. Section fj9] summarizes our main re- 
sults and conclusions. A comparison of our LF measure- 
ments with published observational estimates is shown 
in Appendix [XI 

2. DATA 

The H UDF ACS WFC data are the same as de- 
scribed bv lBeckwith et all ()2006D and include exposures 
in the F435W (hereafter B435 ), F606W (hereafter Veae), 
F775W (hereafter 1775), and F850LP (hereafter zgso) fil- 
ters for a total of 400 orbits. We used the publicly re- 
leased HUDF images without further processing. 

As a part of the UDF05 project, ACS WFC data were 
obtained for the NICP12 field in the F606W, F775W, and 
F850LP filters amounted to a total of lOI HST orbits; Ta- 
ble [1] lists the specifics of these observations. The center 
of the NICPI2 field is located at RA(2000):=03:33:03.60, 
Dec(2000)=-27:41:01.80; Figure [1] shows aU UDF05 
pointings relative to those of HUDF ACS and NICMOS- 
parallels. The UDF05 ACS observations of the NICPI2 
field employ a larger dithering step than those adopted 
for the original HUDF; this choice was made in order to 
obtain an improved sub-pixel sampling for the NICMOS 
parallel data that were acquired simultaneously, relative 
to the original HUDF NICMOS-parallel observations. 

The exposure times for the UDF05 observations were 
designed so as to achieve a constant AB-magnitude depth 
in all filters for all fields. The individual ~1200s expo- 
sures (two in each orbit) were combined and re-pixeled to 
a scale of 30 mas using the task " multidrizzle" ; the task 
"tweakshifts" was used to optimize the shifts between dif- 
ferent images. Supcrdarks and superbiases were obtained 
combining several hundred single frames to increase the 
S/N. Finally, a herringbone artifact was removed from 
the data by processing the images with a Fourier filter, 
and the electronic ghost present in the ACS images was 
eliminated by means of an empirical model. The model 
is built upon the fact that the electronic ghost of a source 
is located in quadrants different from the one containing 

^'^ Throughout this paper we adopt the concordance cosmology: 
flM = 0.3, = 0.7, Ho = 70 kms -^Mpc-\ i.e h = 0.7. Magni- 
tudes are given in the AB system llOke fc GunnllTOSg) . 
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TABLE 1 

Properties of the UDF05 and HUDF Observations 



Field 


Filter 


HST Orbits 


Exposure Time [s] 


lOcr-magnitude^ [AB] 


NICP12 


Vsoe 


9 


21'600 


28.49 


NICP12 


«775 


23 


54'000 


28.44 


NICP12 


2850 


69'' 


168'000 


28.47 


HUDF 


■B435 


56 


134'900 


29.38 


HUDF 


Veo6 


56 


135'300 


29.81 


HUDF 




144 


347'100 


29.43 


HUDF 


2850 


144 


346'600 


28.73 



^ within an aperture of O'.'IS radius'' 70 planned; one lost due to loss of lock 



the source and on the fact that ACS/WFC is read out 
from each corner so that if one dithers a source closer 
to its readout corner the ghost will move in the opposite 
direction closer to its own corner. Thus, when we drizzle 
together a dithered pair of images, we align the sources, 
and smooth out the ghosts. Our algorithm (conceived by 
L.E. Bergeron) is based on the idea of flipping each quad- 
rant before drizzling. This simple step coadds the ghosts 
and smooths out the sources and, when applied itera- 
tively with proper masking, can be used to obtain clean 
images of the ghosts which are then subtracted from the 
images. These last two steps were developed after the 
first public release of the UDF05 data (version 1, "vl"). 
The newly processed images of the UDF05 that we use 
in our analysis represent an improved "version 2" (v2) of 
the final products. 

The PSF FWHM of the UDF05 final ACS images are 
about O'.'l in aU passbands. In the UDF05 NICP12 im- 
ages, the 10 a limiting magnitudes in apertures of 0'.'15 
radius are 28.49, 28.44, and 28.47 respectively in Vqqq, 
1775, and Z850 (see Tabled]), i.e. the planned uniform 
sensitivity in all passbands is achieved. For the HUDF 
the i?435, Vgoe, 1775, and zgso 10 c limiting magnitudes 
in apertures of 0'.'15 radius are 29.38, 29.81, 29.43, 28.73 
respectively. 

2.1. Noise Correlation 

A point kernel was used for all passbands when driz- 
zling the data, with the exception of the NICP12 Vgoe 
image, for which the pixels were not shrunken to a point 
before drizzling to the output frame (pixfrac ^0), as the 
dithering pattern was inadequate to properly sample the 
pixels. As a result, this image is expected to show signifi- 
cant noise correlation. However, the noise level measured 
by computing the rms of cleaned, block averaged areas 
of the Vgoe images, or by computing the total flux in the 
auto-correlation peak shows a higher degree of correla- 
tion than the one due to ditherin g which is well ch a racter - 
ized by the model published bv lCasertano et al] (|2000D . 
Moreover, some noise correlation is also measured in the 
^775 and ^850 images which were not expected to have any. 
The unexpected correlated noise seen in these images is 
of the same amplitude as the additional component seen 
in the Veoe image and we believe that this residual noise 
correlation is introduced by the reference files (flat fields, 
biases, darks) and, possibly, by unresolved background 
sources. Any feature present in a given pixel in the ref- 
erence files is propagated to a larger area on the final 
science images by the adopted dithering pattern. Ne- 
glecting this noise component would lead to overestimat- 



ing the depth achieved by the observations; for example, 
in the NICP12 field, the noise would be underestimated 
by 30% for galaxies with an area of 0.3 x 0.3 arcsec^. 

Relative to our final ACS UDF05 images, the released 
ACS HUDF images, used in this paper, are more severely 
affected by correlated noise, due to presence of the un- 
corrected herringbone noise and electronic ghost in ACS 
(as these effects had not been characterized at the time) . 
This explains why, despite the factor of about two in 
exposure time, the HUDF zgso image is only 0.26 mag 
deeper than our NICP12 zgso image. 

The rms maps produced by multidrizzle were rescaled 
in order to match the expected flux errors of galaxies with 
a measured area of 80 pixels, which corresponds to the 
median size of our dropout galaxies. No simple renormal- 
ization is optimal for all source sizes. The adopted choice 
underestimates the S/N for point sources and overesti- 
mates it for sources that are more extended than the area 
considered; however, it provides on average, an optimal 
renormalization for the typical galaxies we study in this 
paper. Simulations are used to assess proper errors as a 
function of the sources' physical parameters. 

2.2. Source Detection 

Source detection and photometry measurements 
were done with the software tool SExtractor 
(jBertin fc Arnouti Il996( l. which was run in double 
image mode with zgso as the " detection band" . Despite 
the fact that the 1775 images for both the HUDF and the 
NICP12 data reached a similar, if not deeper, magnitude 
limit than the 2350 band, the latter was preferred for 
the source detection because of the larger uncertainties 
in an i775-selected LF, which are driven by redshifting 
of the Lya forest into the 1775 passband at z > 4.8, 
significantly decreasing the observed fluxes. 

The SExtractor parameters were optimized so as to 
maximize the number of detected galaxies while mini- 
mizing the number of spurious sources. For NICP12, the 
detection threshold was set to be a minimum of 10 con- 
nected pixels 0.55 a above the background. This led 
to 5,734 objects within the central, highest S/N 11.9 
arcmin^ of the zgso NICP12 image. Similar SExtractor 
configurations returned 6,336 objects in the central 11.2 
arcmin'^ of the HUDF 2:350 image. 

In order to estimate the contamination of spurious 
detections in the source catalogs, we cleaned the im- 
ages from the detected sources, and reran SExtractor us- 
ing the same configuration parameters on the " inverted" 
(i.e., multiplied by -1) images. This contamination was 
found to be < 0.5%. A final visual inspection was per- 
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observed wavelength [Angstrom] 

Fig. 2. — Synthetic SED of a starbursting galaxy at z = 5 
with ^850 = 25.5 mag. The optical passbands used for the UDF05 
program are overplotted. The vertical dashed lines indicate the 
hya edge at redshifts 4.5, 5 and 5.5, respectively. 

formed to clean the catalogs of residual spurious sources 
and stellar diffraction spikes. The source magnitudes 
were corrected for the small amount of Galactic extinc- 
tion expected in the field, E(B-V ) = 0.007 an d 0.008 in 
NICP12 and HUDF, respectively (jSchlegel et a l. 1998). 

Colour measurements were derived from the galaxy 
ISO magnitudes (apertures matched to the detection 
isophotes), since these are closely matched to the higher 
S/N parts of the objects and, as shown by our simulations 
(Section 2]), include much less noise than the SExtractor 
AUTO magnitudes which we use for the total magnitude 
of a source in the detection band. This is especially im- 
portant for the very compact and faint sources that we 
study in this paper. 

3. THE SELECTION OF Z - 5 LYMAN-BREAK 
GALAXIES 

3.1. The Color- Color Selection Criterion 

The redshift evolution of synthetic galaxy models, 
properly attenuated by intergalactic hydrogen absorp- 
tion, can be used to identify regions in a color-color di- 
agram which efficiently disentangle galaxies in a specific 
high redshift window from lower-redshift interlopers. 

For the identification of a Vgoe-dropout sample (see 
Figure [2l) , this selection window has been optimized to 
be (see iGiavalisco et al... 2004b: Bcckwith et al.ll2006l ): 

Veoe - i775 > min [2, 1.5 -I- 0.9 • (1775 - zgso)] (1) 
^606 - i775 > 1.2 (2) 

«775 - 2^850 < 1-3. (3) 

Galaxies in the range z ~ 4.5 — 5.7 are well identified by 
this color-color criterion; the reliability of the method has 
been shown to be very efficient (^ 90%) by spectroscopic 
follow-up surveys (|Vanzella et al.|[2006f ). 

The elimination of low redshift interlopers is facilitated 
by the availability of the -B435 passband for the HUDF 
sample. This is however not available for the NICP12 
data. We therefore used the HUDF data to check the im- 
pact of including or excluding t he additional S/N ( i34.ci5)< 
3 cr iterion that was adopt ed by IGiavalisco et ahl (|2004bD 
and lBeckwith et all ()2006f ) in order to minimize the con- 
tamination from low-z interlopers. Only one object was 
rejected from our HUDF dropouts catalog on the basis of 
this _B435-band constraint. Therefore, for uniformity with 
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Fig. 3. — Evolutionary tracks of different types of galaxies, and 
location of the stellar sequence, in the Veo6 ~ *750 versus 1775 — zsso 
color-color diagram. The tracks are constructed using the Bruzual 
& Chariot (2003) population synthesis models. The thick black line 
corresponds to the adopted selection (Eq. [T]-|3]l. The thin solid 
lines show the tracks of 200 Myr old, continuously star-forming 
galaxies at 2 > 4 with E ( B — V ) = 0,0.15,0.3 (indicated in the 
figure) for a lCalzetti et al.l I I2000I ) dust extinction relation. The red- 
shifts steps 4.5, 5, and 5.5 are indicated with small dots. Lower red- 
shift galaxy types are plotted with dotted, dash-dotted and dashed 
lines [z = — 4), corresponding to local irregulars, Sbc spirals, and 
a single stellar population with an age of 3.5 Gyr, respectively. The 
latter was chosen because it corresponds to a maximally old stellar 
population at 2 ^ 1.5, the redshift at which single-burst early-type 
galaxies contaminate the selection . The s mall crosses correspond 
to the 131 Galactic stars from the IPickle j ||T998[) library. 



the analysis of the NICP12 field, for which the is not 
available, we did not include any i?435-band constraint in 
our selection of the HU DF Vaoe -dropo uts sample. 

The stellar library of iPicklesI ()1998 f) shows that stars 
are expected to lie within a very well defined region on 
the Vgoe — *775, *775 ~ ^850 diagram. Indeed, the few ob- 
vious stars in our field lie within this area. To further 
optimize the removal of stars from our catalogs, the SEx- 
tractor stellarity parameter was also considered. This 
parameter is determined by application of a neural net- 
work approach to estimate the likelihood that an object 
is a point source. However, this method is not very reli- 
able at faint magnitudes, and therefore cannot be used to 
identify stars close to the magnitude limits. Inspection 
of the Vgoe ~ *775 vs. 1775 — zgso diagram and of the half- 
light radius (fi/2) vs. zgso diagram showed that sources 
with stellarity > 0.85, zsso < 27 and ri/2 < O'.'l are most 
likely stars in the NICP12; hence all such objects were 
excluded from our analysis. A similar criterion was ap- 
phed to the HUDF for zgso < 28. Only two and one of 
such sources were within our Vgoe-dropout selection win- 
dow for the NICP12 and HUDF samples^^, respectively; 
these objects however had already been visually identi- 
fied as stars. Although unlikely, stars fainter than 2359 
~ 27 and 28 may remain a (small) source of contamina- 
tion in our NICP12 and HUDF catalogs, respectively. 

In Fig. [3] the evolutionary tracks of different galaxy 
types are shown in the Vgoe — ^750 versus 1775 — zs5o color- 
color diagram. The tracks are built using the population 
synthesis models of Bruzual & Chariot (2003). The at- 
tenuation by intergalactic hydrogen causes essentially all 

A ll other sources in the HUDF stellar catalog of lPirzkal et al.l 
II2005I ) that we do not identify as stars lie outside and very far from 
the border of our color-color selection window. 
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the galaxy types at redshifts above ~4.5 to enter the se- 
lection window. Elliptical galaxies at z ^ 1.5 are a source 
of contamination in the selection criterion due to their 
red colors; this contamination is however expected to be 
small in our Vgoe-dropout samples, given the expected 
low surface densities of these galaxies. In fact, none of 
the seven spectroscopically identifie d passively evolving 
galaxies at z > 1.4 in the HUDF from lDaddi et all ()2006h 
were selected as Vgoe-dropouts. 

3.2. Contamination from Low-Redshift Interlopers 

We have estimated the possible fraction of interlop- 
ers by applying our selection criterion of Eqs. (1)- 
(3) to a library of ^ 3000 synthetic SEDs built upon 
Bruzual-Charlot (2003) models, adopting the LF de- 
rived bv Steidel e t al. (199ft ) at z ~ 3 and no evolution. 
The models includ e the effects of intergalactic absorp- 
tion (jMadaulll995l ). and span a wide range of metalhci- 
ties (0.04 — 2. 5^0), dust reddening, emission lines, and 
different star formation histories. These included simple 
stellar populations, models with continuous star forma- 
tion at constant metallicity, models with self-consistent 
enrichment during the continuous star formation (esti- 
mated with a closed box model and with a model with 
infall of zero- metallicity gas), and two-burst models com- 
bining with five different mass ratios an old (> 50 Myr) 
stellar population with a young (1 — 10 Myr) star form- 
ing population that includes self-consistent nebular line 
and continuum emission. For each model, versions with 
and without dust reddening were computed; eight val- 
ues of reddening were implemented assuming a galac- 
tic dust screen with Ay logarithmically spaced between 
0.05 and 6.4. Emission lines and nebular continuum were 
computed in the approximation of case-B recombination 
for Hydrogen and Helium, based on the ionizing flux 
from the stellar SED; Oxygen and Carbon line inten- 
sities were derived from an analytical interpolation of a 
grid of calculations performed with the code CLOUDY 
(|Ferland et al.|[T998[) . Gas temperatures were derived as 
a function of metallicity. 

The resulting redshift distribution is shown in Fig. [4l 
The non-uniform distribution of the histogram is due to 
SEDs that satisfy the selection criterion only for por- 
tions of the redshift interval. For instance a 1 Gyr old 
solar metallicity model with modest extinction satisfies 
the criterion only for z > 5.5. We identify three possible 
classes of low redshift interlopers at z ~ 0.6, ~ 1.6, and 
~ 3.2. They are generally very reddened SEDs with ei- 
ther emission lines or an underlying old population. The 
total interloper fraction from these objects is found to be 
23% (17% if one considered only z < 4 as interlopers). 
This is only a rough estimate, as it does not account 
for the different galaxy populations having different and 
possibly evolving LFs. However, previous works show 
that the utilization of this library overpredicts th e num- 
ber o f interlopers (e.g. with respect to Malhotr a et al.l 
|2005() : also, in the spectroscopic follow-up survey of the 
Vgoe-dropout galaxies in t he Great Observatorie s Ori- 
gin Deep Survey (GOODS: iGiavalisco et aHl2004a^ ■ only 
^ 10% of stars and l ow-z galaxy interlopers were found 
(jVanzella et al.ll2006( l. Therefore, the real fraction of in- 
terlopers is possibly smaller than our conservative esti- 
mate. Note however, that we do not account for this 
source of systematic error in the rest of our analysis. 
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Fig. 4. — The predicted redshift distribution for Veoe-dropout 
-dropouts as derived assuming a non-evolving LF and synthetic 
SEDs (see text). The total interloper fraction is estimated to be 
23% and is primarily contributed by very reddened galaxies with 
either strong emission lines or an old stellar population component. 
The inset panel zooms over the primary redshift range of interest. 
The non-uniform distribution of the histogram highlights the fact 
that different SEDs are responsible for the signal in different red- 
shift intervals. 

4. CATALOGS COMPLETENESS AND 
SIMULATIONS 

The catalogs produced by SExtractor are subject to 
certain biases, which depend strongly on the chosen pa- 
rameters. In addition to the detection settings, i.e., num- 
ber of connected pixels and sigma-threshold above the 
background, the parameters that control the deblending 
and cleaning of sources are of particular importance. The 
correct interpretation of the LP that is obtained from the 
SExtractor-based catalogs requires a quantitative assess- 
ment of these biases. To this end, we performed two in- 
dependent sets of simulations, the first based on idealized 
galaxy profiles and the second on the observed galaxies 
themselves. The simulations are described in detail in 
Appendix [Bl Briefly: 

1. In the first set of simulations (which we refer 
to as the " idealized-galaxies" simulations), we 
constructed an extensive set of exponential- and 
deVaucouleurs-profile galaxies, that were appropri- 
ately PSP-convolved and inserted at random posi- 
tions into the NICP12 and HUDF images. SEx- 
tractor was then re-run with identical parameters 
as for the compilation of the original data catalogs, 
and the resulting object catalogs were used to iden- 
tify detected, undetected or blended sources. This 
approach is similar in spirit to what was done by 
the GOODS team (|Giavahsco et al.ll2004al ). 

2. In the second set of simulations (which we refer 
to as the " dimmed-galaxies" simulations), we gen- 
erated a set of test images by co-adding to each 
original image a copy of itself, shifted by a few pix- 
els and dimmed by a given amount. The composite 
images were rescaled to have the S/N vs. magni- 
tude relation as the original image, and SExtrac- 
tor was run on these final test images using the 
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same parameters as for the extraction of the orig- 
inal data catalogs. This second set of simulations 
provides a complementary analysis to the above, as 
it offers the advantage of working with real rather 
than idealized galaxies, and does not depend on 
uncertainties in the adopted PSF. 

4.1. Input versus Output Properties 

A problem in establishing the completeness of the 
source catalogs is that the true magnitudes of the sources 
are expected to be typically brighter than the measured 
ones, due to losses of flux from the wing s of the galaxy 
light distribution. ISirianni et al.l HooSy measured the 
flux losses for point-sources in the ACS camera fllter- 
bands and found that they strongly depend on the SED 
of the source, due to a stronger scattering halo of the 
PSF at longer wavelengths. Those authors find offsets of 
~ 0.1 mag for Veoe and 1775, respectively, and ~ 0.15 mag 
for Z850 in apertures of radius 0.4". To take into account 
the extended nature of our catalog objects, we used our 
two sets of simulations to characterize the relation be- 
tween input (true) parameters and output (measured) 
quantities. 

Idealized deVaucouleurs- and exponential-profile 
galaxies were generated with random zs5o magnitudes 
between 24 to 31. Two separate sets of tests were 
performed, the first using sizes^^ uniformly distributed 
in the range 0.05 to 0.5 arcseconds, and the second using 
sizes distributed according to a lognormal distribution 
centered at 0.25 arcseconds with a = 0.3. The results 
were essentially independent of the adopted input size 
distribution. In the following, we only use the simula- 
tions with uniformly distributed input sizes. Measuring 
the source magnitudes within an elliptical aperture of 2.5 
Kron radii (SExtractor AUTO mag), compact sources 
were detected with 50% probability down to 28.75 and 
29 magnitudes in the NICP12 and in the HUDF images, 
respectively. As expected, the detected fluxes were 
typically smaller than the input ones. By integrating 
the theoretical deVaucouleurs and exponential proflles 
out to 2.5 times their Kron radii, the theoretically 
expected flux losses are found to be 9.6% and 4.0%, 
respectively, corresponding to 0.11 and 0.04 magnitude 
offsets. However, since the Kron radii measured by 
SExtractor are also underestimated, the offsets are 
expected to be even larger than those above. Indeed, in 
our simulations we found a magnitude-dependence for 
the amount of dimming, from ~ 0.3 mag and ~ 0.07 mag 
at 2:350 — 24.25 to ^ 0.6 and ~ 0.2 mag at zgso = 28.25 
for the two kinds of proflles respectively. 

We also used the dimmed-galaxies simulations to ob- 
tain an independent estimate. In this case, the differ- 
ences between dimmed and original magnitudes are ex- 
pected to be smaller than in the case of idealized galaxies. 
The offset was indeed marginally smaller than for the 
exponential-profile idealized galaxies, i.e., in the range 
Amag ~ 0.04-0.12. 

For very compact galaxies {ri/2 £ 0'.'15), the measured 
sizes were generally overpredicted due to PSF blurring, 
while the opposite was true for larger objects. These 
biases showed a strong dependence on magnitude. The 
faintest idealized deVaucouleurs galaxies were typically 

In all our analysis size always refers to half-light radius. 



measured to have ^ O'-'l (i-e. to be unresolved at the 
resolution of the HST), independent of their original size. 
This offers a challenge for studying the size evolution of 
LBGs with redshift. 

Furthermore, since we measure the colors by relying on 
the SExtractor dual-image mode, i.e. using the zsso de- 
tection apertures, it is not adequate to apply a constant 
offset to the Vgoe and 1775 magnitudes. The color mea- 
surements depend in fact on the zgso flux as well as on 
the intrinsic color. Both the 1775 — zs5o and the Vgoe^iyrs 
colors are affected in such a way that red sources have 
even redder, and blue sources even bluer, measured col- 
ors. This is not including the dependence of the PSF 
halo on the color. 

In the light of these complex and uncertain corrections 
for dimming effects, we chose not to implement any cor- 
rection to the measured catalog magnitudes. Thus, we 
note that our magnitude measurements might be under- 
estimated by an amount, as estimated from our tests, of 
up to about 0.2 magnitudes. We found, however, that 
these light losses had only a marginal effect on our final 
LF when extending it with the SDF (see section [5^ : the 
resulting slope is steepened by about Aa ~ 0.03. 

4.2. Detection Completeness 

There are two possible reasons why galaxies could be 
missing from our catalogs. First, a galaxy might not be 
detected due to its too low surface brightness; second, it 
could get blended together with another source and pro- 
duce a false catalog entry. These two effects are strongly 
dependent on the SExtractor parameters. 

We used both the idealized-galaxies and dimmed- 
galaxies simulations to estimate the strength of these ef- 
fects in our catalogs and final samples of Vgoe-dropout 
galaxies. Specifically: (a) We associated idealized- 
galaxies with similar output (SExtractor) magnitudes 
and sizes to real galaxies randomly selected from our 
original catalog. We then estimated the completeness 
corrections for each magnitude bin by comparing the 
number of sources in each bin of SExtractor magnitude 
with the total number of idealized-galaxies of different 
input (i.e., theoretical) magnitudes and sizes that con- 
tributed to that specific bin of output magnitude, (b) 
Using the dimmed-galaxies simulations, we computed 
the magnitude-dependent completeness corrections by 
re-scaling the number of galaxies in each bin of output 
magnitude (i.e., recovered for the dimmed galaxies by 
SExtractor) to the number of sources of different input 
magnitudes (i.e., the theoretical values of the dimmed 
magnitudes) that contributed to that specific bin of out- 
put magnitude. Figure [5] shows the results of our tests. 
For both fields, the two independent estimates for the 
completeness in our samples, derived from the idealized- 
and dimmed-galaxies simulations, are in good agreement. 
At the magnitude limits of our selections, the NICP12 
and HUDF catalogs are expected to be about 65% and 
50% complete, respectively. 

The probability that a given source is blended by SEx- 
tractor with another object depends on the magnitude of 
the source. In our simulation tests, when the SExtractor 
detection blended a simulated source with a real galaxy, 
we considered the simulated object as a detection when 
this was brighter than the original galaxy. Figure [5] also 
includes the incompleteness due to source blending; this 
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TABLE 2 

Number of detected LBGs and corrections in the 
two fields 



mSSo'' Ncat^ Ncorr" Cldcalizcd'^ CDimmcd'' 




Fig. 5. — Estimated detection completeness over the range of in- 
terest in the two fields as a function of observed magnitude from the 
idealized-galaxies (solid line) and dimmod-galaxies (dashed line) 
simulations. The completeness docs not reach 100% even at the 
brightest magnitudes, which is due to blending of sources by the 
SExtractor algorithm (dotted line, for idealized-galaxies). 

is the reason why even our brightest sources are only 
complete at the 95% level. 

5. THE V^ioe-dropout CANDIDATES 

The application of the selection color window discussed 
above to our NICP12 and HUDF catalogs down to 28.5 
and 29.25 magnitudes (with a S/N > 6 cut) produces 
a total of 101 and 133 candidate z ~ 5 LBGs, respec- 
tively. Note that for sources whose Vgoe ISO fluxes were 
fainter than their 2(T-errors, we replaced the measured 
Vgoe fluxes with the corresponding 2a values in order to 
get a reliable lower limit on their Vgoe — *775 color. 

To correct our catalog for both random and systematic 
photometric errors we adopted a Monte-Carlo technique. 
The amplitudes of the errors were estimated from both 
sets of simulations described in section 31 The simula- 
tions show that systematic color shifts (i.e. input ver- 
sus recovered color) are independent of the galaxy pro- 
file, unlike offsets in individual passband magnitudes. 
Therefore, before the application of random errors, the 
Vgoe— «775and 1775 — zssocolors of each catalog source were 
corrected for systematic errors; only a random error was 
instead applied to the zs5o magnitude (see also section 
14. ip . The selection of Vgoe-dropouts was iteratively re- 
peated and, after of order several hundred iterations, the 
mean number of Vgoe-dropout sources was computed in 
bins of Z850 magnitude. Table[2]lists the mean number of 
Vgoe-dropout sources Ncorr per magnitude bin for both 
our fields which we use in the remainder of the analysis. 

Tests were performed to ensure that the final LBG LF 
was independent of the adoption of the 2a Vgoe flux limits 
for faint Vgoe sources. In particular, the calculations were 
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^ central bin magnitude'' uncorrected number of color- 
selected LBG sources'^ corrected number of Veoe-dropout 
sources after Monte-Carlo resampling'^ completeness de- 
rived from the idealized galaxies simulations'^ complete- 
ness derived from the dimming simulations 



repeated adopting la and 3a limits for the flux substi- 
tutions; no significant differences were found in the final 
LFs, as the probability of selecting a galaxy at any given 
redshift changed accordingly (see next section). 

6. THE LUMINOSITY FUNCTION OF 
LYMAN-BREAK GALAXIES AT Z - 5 

The rest-frame UV-continuum LF was derived from a 
maximum likelihood fit of a Schcchtcr (1976) function to 
the observed number of Veoe-dropout galaxies apparent 
magnitude bins, Nf. 

Ni= dm N{m;(j)^,M^,a) ^ (4) 

J mi 

f°° dV /■"" 

/ dz—— / dm p{m,z)(j){M{m,z);(j)^,M,,,a) 

JQ "2 Jmi 

Here mi and m„ are the lower and upper bounds of 
the z-th magnitude bin, centred around m.^. The factor 
p(m, z) is the probability that a galaxy with magnitude 
m at redshift z is detected and selected to be a 2: ~ 5 
LBG. 

An "effective volume" (Vcs) technique was also used to 
compute the stepwise LF in bins of absolute magnitude 
under the assumption that M{m, z) is slowly-varying 
with redshift, i.e.: 

N,~(t>{M{m,,z)) / dz—p{m„z) (5) 

= 0(M(to„z)) Ueff(m,) (6) 

Over the redshift range of interest, z — 4.5—5.7, the abso- 
lute magnitude at any given apparent brightness changes 
by about 0.5 magnitudes, i.e., an amount comparable 
with the size of our magnitude bins. The resulting LF is 
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thus generally less robust than the one computed using 
Eq. m 

In order to properly constrain the value of in our 
best fit Schechter function, we adopt this v e rsion of our 
LF for matching it to the lYoshida et al.l (j2006) data 
which constrained the z ~ 5 LF at bright magnitudes 
in the Subaru Deep Field (SDF) (see Sectioning. 

The K-corrections required to transform the observed 
magnitudes into rest-frame 1400 A absolute magni- 
tudes, M{m, z), were computed assuming a 200 Myr old, 
continuously starformi ng galaxy template, reddened by 
E{B-V) = 0.15 with a lCalzetti et"air(|200(]( ) dust extinc- 
tion relation (hereafter referred to as "starburst" dust ex- 
tinction). The differences in the K-corrections for alter- 
native reasonable assumptions on the stellar population 
properties, including significantly younger stellar popu- 
lations, were checked and found to be negligible (< 0.05 
mag). 

The probability that a galaxy with magnitude m at 
redshift z is detected and selected to be a LBG at z ~ 5 
can be decomposed into the product of two compo- 
nents, i.e., the completeness-probability C{m) and the 
selection-probability S(m,z), i.e. p{m,z) = C{m) ■ 
S{m, z). We discuss this below. 

6.1. The Selection Probability S{m,z) 

The selection-probability as a function of magnitude 
and redshift was computed including the scatter in the 
photometry, which was derived from our simulations, and 
replacing with the 2 — a flux value the Vgoe magnitudes 
fainter than this limit. The input assumptions for the dis- 
tribution of SEDs in th e Vfine-dropout popu lation were 
based on the findings of ISteidel et al.l ()1999| ). These au- 
thors inferred a distribution of spectral slopes from their 
sample of spectroscopically-confirmed z ~ 3 LBGs by 
assuming a 1 Gyr old, continuously starforming galaxy 
SED with varying dust content acc ording to the starburst 
dust extinction relation (see also lAdelberger fc Steidel 
I2000f) . The E{B - V) distribution was found to be 
roughly approximated by a Gaussian with a median of 
0.13 and a cr = 0.12. Note that this results in a small 
fraction of the galaxy population which is fitted by non- 
physical, negative values of E{B—V)] this might indicate 
the presence of strong Lya emission for some galaxies. 
Nevertheless, this "effective" distribution of E{B — V) is 
in good agre ement with the typical dust attenuation val- 
ues found bv iPapovich et al] (|2001h . who used synthetic 
SED fitting to z ~ 3 LBGs to infer the properties of these 
galaxies. Since an age of 1 Gyr might be too old for z ~ 5 
galaxies, we used similar models with a 200 Myr old , con- 
stantly starforming SED (jBruzual fc GharlotH'2003D . The 
selection probability was found to be virtually indepen- 
dent of the assumed stellar age; this is understandable in 
the light of the dominant contribution from the youngest 
stellar populations to the UV emission. 

To explore the effects of different assumptions for the 
intergalactic hydrogen absorption, the SE Ds were atten- 
uated according to the calculations of both lMadaul ()1995D 

Note that the situation would be exacerbated for an i-n^- 
selected (rather than 2350 -selected) Vgoe-dropout sample, since 
the differences in absolute magnitude over the same redshift range 
could be up to ~2 magnitudes, i.e., much larger than the size of 
the bins, due to intergalactic Lya absorption. 
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Fig. 6.— Top - The NICP12 selection probabilities for the 
different assumptions on the SED distributions and intergalactic 
hydrogen absorption discussed in the text, for the case of a bright 
(25 — 25.5, blue lines) and faint (28 — 28.5, red lines) magnitude 
bin. The solid (dot-dashed) and dashed (dotted) lines refer to the 
starburst (SMC) dust extinction relation, assuming the intergalac- 
tic hydrogen absorption prescription of Madau (1995) and Meiksin 
(2006), respectively Bottom - The same for the HUDF. 

and'Meiksin' (2006). In the latter case, the selection prob- 
abilities are slightly reduced, as a consequence of the 
smaller attenuation and thus smaller reddening of the 
SEDs. 

We estimated the impact on the selection probability of 
assuming that there are no Lyman-Limit systems along 
the line of sight attenuating the observed Vgoe flux, a sit- 
uation which is expecte d to happen w ith a probability 
of only 6% up to z ~ 5 ()Meiksinl 12001 . This effect was 
found to be negligible, i.e., the selection probability was 
reduced by less than 5% around z ^ 5.2. In the remain- 
der of our analysis we therefore assume the estimated 
mean values for the intergalactic absorption. 

We also tested al t ernat ive distributions for the UV 
slopes. IVerma et al.l (|2007D used stellar population syn- 
thesis models to derive ages, masses, and dust atten- 
uation properties of z ^ 5 LBGs. Using the Small 
Magellanic Cloud (S MC) dust extinction l aw, both 
IPapovich et"al1 (|2001l ) and lVerma et al.l (|2007f ) find that 
the typical UV attenuation is smaller than in the case 
of the starburst extinction law. The exact values some- 
what depend on the metallicities and Initial Mass Func- 
tions (IMFs). To estimate the selection probability, we 
adopted a 100 Myr old, continuously starforming tem- 
plate reddened by a Gaussian distribution of E{B — V) 
with mean= 0.05 and a = 0.04, which is approximately 
in agreement with the distributions found by Papovich 
et al. and Verma and et al. Also in this case we explored 
the effects of changing the prescriptions for the hydrogen 
intergalactic absorption. Differences of order 10% or less 
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were found with respect to the previous estimates. 

The selection probabihties for the NICP12 field, as de- 
rived under the various different assumptions, are shown 
in Fig. [6] for a bright and a faint magnitude bin. In gen- 
eral, sinc e the intergalactic absorption relation derived by 
iMeiksinI ((2006ll is smaller than the corresponding calcula- 
tion of Madau (1995), galaxies are less reddened assum- 
ing the former prescription; consequently, at any redshift 
of interest, they lie closer to the boundaries of our color 
selection. This leads to a reduced probability that they 
are selected as LBG candidates. The effect is exacerbated 
when the SMC extinction law is assumed: a consider- 
able drop in the redshift distribution is expected around 
z ~ 5.3, and the lower limit for the selection-redshift is 
slightly increased. Note that the adoption of the 2 — a 
lower limit for the Vgoe fluxes leads to a considerably 
different redshift window for the fainter magnitude bin: 
the reddest, highest-redshift objects are not detected. 
Furthermore, due to the larger photometric scatter at 
fainter magnitudes, the redshift selection smears out to 
substantially lower redshifts, down to z ~ 4. Note that 
we have not applied any correction for incompleteness in 
these plots; therefore, the drop in selection probability 
at fainter magnitudes that is observable in the Figure is 
exclusively due to the errors in the color measurements. 
Similar results hold for the HUDF, for which however 
the selection probabilities are somewhat larger, due to 
the smaller uncertainties in the color measurements (see 
bottom panel of Fig. [6|). 

6.2. The Impact of Different Physical Assumptions on 

the LF 

As we discuss below, different assumptions for physical 
input parameters such as the dust extinction law and the 
hydrogen intergalactic absorption relation have a non- 
negligible impact on key quantities such as the effective 
volume and thus the average luminosity density. It is 
thus of paramount importance that LFs derived from dif- 
ferent data sets are first rescaled to identical assumptions 
for SED distributions, dust properties and intergalactic 
extinction, before they can be compared to study any 
possible evolutionary effect in the galaxy populations. In 
our analysis we adopt the Madau (1995) and starburst 
relations as our "fiducial" assumptions for the hydrogen 
absorption and dust extinction, respectively, as these are 
widely used in most previous works and allow for a direct 
comparison with the published results. We explore how- 
ever below the effects on the LF parameters of varying 
this assumptions. 

Our pencil-beam observations do not properly con- 
strain the value of in the LF: given the rarity of 
very bright objects, the small volumes probed by the 
NIC? 12 field and the HUDF are not sufficient to ade- 
quately sample the knee of the Schechter function. We 
therefore rely on the SD F z ^ 5 LF of Yoshida et al. 
(2006; see also Appendix IA.1[) to constrain the bright- 
end of our Vcs LF. These authors used similar fiducial 
assumptions for the intergalactic absorption and dust ex- 
tinction relations, and similar techniques to correct for 
the observational biases, to those adopted in this paper. 

The SDF LF is based on the Ri'z' SUBARU filters, and 
is sensitive to a redshift window that is only slightly nar- 
rower than that of our data. It reaches down to z' = 26 
mag, i.e., ~ 3 magnitudes brighter than our nieasure- 




FlG. 7.— Total Vcff LF obtained combin ing the NICP12 and 
HUDF data (this work) with the SDF LF IIYoshida et al.l l200l . 
ffi'z'-sample). Red squares and blue circles correspond to the 
NICP12 a nd HUDF V^ff LFs r espectively. Green crosses are the 
Voft-LF bv l Yoshida et an ( |2006i ). The errorbars indicate Poissonian 
errors only. The solid line shows the best-fit Schechter function to 
the total LF (see text and Table [Sj. 

ments. In combining the SDF and NICP12/HUDF LFs, 
we neglected small K-corrections from 1500 A to 1400 A; 
the Voff SDF LF nicely matches our Kff NICP12 LF in 
the region of overlap. The total Vcs LF, that merges the 
SDF estimate with our NIC? 12 and HUDF LFs, spans 
over the 5-magnitudes range A/1400 = [—22.2,-17.1] 
(Figure ED. 

Adopting the completeness corrections derived from 
the idealized-galaxies simulations, the lowest x^, best 
fit of a Schechter function to the total Vcs LF gives: 



0.9lo 3 ■ 10"^ 
4+0.10 



Mpc -^mag ^, 



-20.78 



+0.16 
-0.16 



mag, a = -1.54:^o;|g (see Figure [7] and Table [3]). As ex- 
pected, the value of is largely influenced by the SDF 
LF, while the slope a is constrained by our fainter data. 
The adoption of the completeness corrections derived 
from the dimmed-galaxies simulations, and of different 
SED distributions, dust extinction laws and intergalac- 
tic absorption prescriptions, has no significant impact on 
this result. 

In particular, all tests performed by varying these 
quantities, independent of whether they include one or 
both of our two deep fields, return a value of that, 
within the errors, is consistent with = —20.7. We 
therefore fixed to this value in order to highlight the 
effects of varying the input physical parameters on the 
LF computed according to Eq. (HI for the HUDF and 
NIC? 12 only. 

The results obtained with fixing the value of Af* = 
—20.7 and excluding the SDF from the fits are listed 
in Table [H Only marginal differences are generally ob- 
served. Within the errors, the derived parameters are 
consistent with each other. Note however that the aver- 
age luminosity density increases by about 20-30% when 
the V606-dropout population is described by the SMC 
dust extinction law and the Meiksin intergalactic absorp- 
tion relation, relative to the more widely used combina- 
tion of prescriptions which we have adopted as our fidu- 
cial model. This is due to a change of a similar order of 
magnitude in the resulting effective volumes. We note 
that the uncertainties in the completeness, estimated 
from our two complementary sets of simulations, have 
a comparably small effect on the faint-end slope of the 
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TABLE 3 
2; ~ 5 Luminosity Functions 



Authors 


Field 


Area 
[arcmin^] 


Mim 

[mag] 


[Mpc 3 mag -■-] 




Ah 
[mag] 


a 


Beckwith et aL (2006) 
Yoshida et aL (20061 
Yoshida et aL (2006) 


HUDF 
SDF 
SDF 


11 
875 
875 


-17.2 
-20.3 
-20.3 


1 ■ 10^3 
O.bStolt ■ 10- 
1.23tlit ■ 10- 


-3 
-3 


-20.5 

-2i.o9is;-?i 

-20.72l»:« 


-1.6 (fixed) 
9 .1-1+0.68 
^••^^-0.60 

-1.82 (fixed) 


This work" + Yoshida 
This work'' + Yoshida 


HUDF+NICP12+SDF 
HUDF+NICP12+SDF 






0. 9 ±0.3- 10" 
1.0 ±0.3 -10" 


'3 
'3 


-20.78 ±0.21 
-20.76 ±0.18 


-1.54 ±0.10 
-1.62 ±0.08 



" Adopting the completeness correction derived from the idealized-galaxies simulations, the starburst dust extinction 
relation, and Madau (1995) intergalactic absorption prescription.'' Adopting the completeness correction derived from the 
idealized-galaxies simulations, the SMC dust extinction law, and the Meiksin (2006) intergalactic absorption prescription. 



LF as the different physical assumptions. 

7. THE ROLE OF COSMIC VARIANCE 

7.1. The Underdensity of the HUDF 

On average, independent of the assumed input param- 
eters, the LF derived from the HUDF resuhs in a lower 
average luminosity density: The HUDF appears to be 
"under-dense" with respect to the NICP12 field (see also 
Fig. [7|). By integrating the LF from Table [4] over the 
magnitude range of overlap between the HUDF and the 
NICP12 field, i.e., M1400 = [-21,-18], the amount of 
this underdensity is about 30%, independent of the spe- 
cific input physical assumptions. Although a larger num- 
ber of low-redshift interlopers in the NICP12 sample can- 
not be ruled out, (a) the homogeneity in the procedures 
with which both data sets have been handled and an- 
alyzed, and (b) theoretical expectations concerning the 
amount of scatter amongst volumes of order of those sam- 
pled with the HUDF and NIC? 12 field, both support the 
interpretation that the under-density in z ~ 5 LEG of the 
HUDF has a physical origin. No te that, in their analy- 
sis of the i-dropouts population, iBouwens et al.l ()2006f ) 
find the HUDF under-dense at the bright end, also by 
about - 30%, relative to the GOODS field. Thus, it 
appears that the ~ 30% under-density of the HUDF ex- 
tends at least over the volume encompassing the entire 
z = 4.5 — 6.5 epoch. 

7.2. Quantification of Cosmic Variance Effects in Our 

Study 

In order to estimate the impact of cosmic variance on 
our measurements, we use a 512'^-particles, dark-matter- 
only cosmological simulation in a box of lOQMpc/h side, 
that results in a particle mass of 5 • lO^M^g/Zi. Simple 
prescriptions to populate the dark matter halos with lu- 
minous galaxies, and a realistic beam-tracing algorithm 
were used to study the field-to-field variations in the 
number counts of star forming galaxies at the epoch of 
interest. 

The numerical simulation wa s carried ou t with the 
pubhc PJVI-Tree code Gadgct-2 ( Springel 2005) adopting 
a (Tg = 0.75 value. The initial conditions w ere generated 
using a code based on the Grafic algorithm (|Bertschingeil 
l2001f ) with a A C DM transfer function c omputed accord- 
ing to the fit of lEisenstein fc Hul (|1999( l. and assuming a 
spectral index n ~ 1. Dark matter halos were identified 
in the simulation snaps hots using the HOP halo finder 
(|Eisenstein fc Hutl[l998[ ): the redshift interval 4.5-6 was 



sampled in steps of Az = 0.125. The simulation repro- 
duces a volume about 73 times larger than the effective 
volume probed by each of our two fields for the Vgoe- 
dropout population 5.7 x 5.7 x 420{Mpc/h)^). 

Source catalogs were constructed for each snapshot, 
and analyzed with our beam-tra cer algorithm, which 
is sirn ilar to the one described in iKitzbichler fc White! 
(2006'). We present details of the beam-tracer algorithm 
in Trenti fc Stiavelll (2007 ). In brief, the pencil beam as- 
sociated with the ACS field of view is traced through the 
simulation box across different snapshots of increasing 
redshift. As the uncertainty in the Veoe-dropout redshifts 
is larger than the size of the simulation box, the simu- 
lation box needs to be beam-traced several times. We 
carried out 4000 Monte Carlo realizations of the beam- 
tracing procedure by varying the initial position of the 
beam; the beam was always angled so as to guarantee 
no overlap between different beam-tracing runs. With 
our choice of angles we got a separation between differ- 
ent beam regions of at least ~15 Mpc/h, which corre- 
sponds to ~540 arcsec at z = 5. We verified that fields 
with such separations are essentially statistically inde- 
pendent. Note also that our lOOMpc/h box is different 
from an ideal simulation of a 420Mpc/h box as in the 
power spectrum we are missing contributions from wave- 
lengths greater than lOOMpc/h. For the scales we are 
interested in at z > 4, the missing power is however neg- 
ligible in the two-point correlation function and hence in 
the cosmic variance. 

Dark matter halos within a beam were flagged, and 
a simple model was used to determine the number of 
Vfjoe-dropouts within a flagged halo. In particular, the 
intrinsic number of galaxies was computed using the halo 
occupation model description. Th is gives an average oc- 
cupation number (jWechsler et al.ll200lh : 

(Ngai) - 0{M - M™„)(1 + (Af/Afi)"), (7) 

where 9{x) is the Heaviside step function, Mmin is a 
minimum halo mass threshold. Mi is the typical scale 
where multiple galaxies are present within the same halo, 
and a > is of order unity. If M > Mmin, in addition to 
the one galaxy at the center of a halo, a suitable number 
of companion galaxies is extracted from a Poisson distri- 
bution with mean (M/Mi)". A fraction 77 of the total 
galaxies in each halo is finally identified as Vgoe-dropouts 
, so as to take into account that Vgoe-dropouts may be 
detected for only a limited period during their evolution. 
In addition to this, we tossed 50% of the selected galaxies 
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TABLE 4 

Impact of Physical Parameters on the LF (Eq. [H Constant M* = -20.7)'!' 



Completeness 


Dust Extinction 


Intergalactic Absorption 


0. NICP12 


a NICP12 


4>* HUDF 


a HUDF 


Idealized 


Starburst 


Madau 


8.4- 10-4 


-1.71 


5.8 ■ 10-4 


-1.74 


Idealized 


Starburst 


Meiksin 


9.5 ■ 10-4 


-1.71 


6.4 ■ 10-4 


-1.75 


Idealized 


SMC 


Madau 


9.0 ■ 10-4 


-1.71 


6.0 ■ 10-4 


-1.75 


Idealized 


SMC 


Meiksin 


11.6 ■ 10-4 


-1.66 


7.0 ■ 10-4 


-1.76 


Dimmed 


Starburst 


Madau 


8.9 • 10-4 


-1.68 


6.0 ■ 10-4 


-1.71 


Dimmed 


Starburst 


Meiksin 


10.0 ■ 10-4 


-1.67 


6.5 ■ 10-4 


-1.72 


Dimmed 


SMC 


Madau 


9.4- 10-4 


-1.68 


6.1 ■ 10-4 


-1.72 


Dimmed 


SMC 


Meiksin 


12.1 ■ 10-4 


-1.63 


7.1 ■ 10-4 


-1.73 



^ Typical errorbars for </>* and a are ±2.1 • 10 4 and ±0.16, respectively. The starburst dust model assumes 
E{B -V) = 0.13 ± 0.12; the SMC dust model assumes E(B -V) = 0.05 ± 0.04. 
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Fig. 8. — Probability distribution of the Veoe-dropouts number 
counts in one ACS field (top) and in the two HUDF + NICPIC12 
fields (bottom), accordingly to different halo occupation distri- 
bution models. Black line: one galaxy per halo, r) = 0.5; red 
line: one galaxy per halo rj = 1\ blue line: HOD model with 
Ml = 5 ■ 10^^ Mq/H and ?y = 1. The Poisson distributions are 
overplotted in green for comparison. 



due to detection and selection incompleteness present in 
our survey. 

The results for different choices of the model param- 
eters are reported in Fig. [51 We consider: (i) rj — 1, 
Ml +00 (i.e., one galaxy per halo); (ii) rj = 0.5, 
Ml — > -f-oo (i.e., 50% probability that a DM halo hosts a 
LBG galaxy); (in) -q ^ I, Mi = 5 • 10"Mo//i and a = I. 
The value of Mmin is kept as a free parameter and it is 
adjusted to have the same average number of counts in 
all three cases. Variations of Mmin are limited in the 
range M™„ = 2 - 3.5 • W^^MQ/h. 

The distribution of the number counts probability is 



very similar in the three cases, and the field-to-field stan- 
dard deviation is cr = 25 ± 2%. This is significantly 
larger than the Poisson uncertainty of 10% for a value 
of ~ 100 counts per field, which is about the number 
of Veoe-dropout counts/field that we have detected. De- 
pending on the specific model, there are minor differ- 
ences in the computed standard deviation (Act « ±2%): 
for 77 < 1, a smaller Mmin is required; therefore, the 
variance decreases, as lower-mass halos become less clus- 
tered. As expected, the largest variance results in the 
case of multiple galaxy occupancy of the halos. 

We also computed the relative variations in expected 
counts for two fields which are separated by about « 10 
arcmin, as is the case for the HUDF and the NICP12 
field. We find that the two fields have essentially inde- 
pendent counts, i.e., the variance of the average number 
of detected Vgoe-dropouts decreases by a factor \/2. 

Finally it is interesting to highlight that the narrow- 
beam geometry of our data significantly reduces the cos- 
mic scatter with respect to a similarly-sized, spherical 
volume for a single field, which has a variance in th e 
number counts of « 50% (see also lSomerville et al.l2004l) . 
This is due to the stretched redshift space of the pencil 
beam, that probes a variety of environments by minimiz- 
ing the probability of enclosing only a large void or a rich 
cluster. 

8. DISCUSSION 

We discuss below the implications of our measurement 
of the faint-end of the z ^ 5 LF of star forming galaxies 
for the evolution of galaxies, in the context of current 
models of galaxy formation. 

8.1. Evolution of the Faint-End of the LBGs 
Luminosity Function between z 3 to 6 

There have been several other studies of the LBG 
LF from z = 3-6. While for z - 3 the dif- 
ferent LFs and number counts see m to be in reason- 
able agreement with each other fe.g. iSteid el et aD ll999l : 
Sawicki & Thompson 2006 ; Hildebrandt et al. 2003), at 
z ~ 6 large ly different LFs we r e reported in the lit- 
erature (e.g. [Bu nker e t all 120041; iDickinson et al.l 120041 : 
lYan fc WindhOT st 200i iMalhotra et"al] 120051 . and see 
also discussion in Bouwens et al. 2005). This might be 
mainly due to the highly model dependent corrections 
needed for the computation of the i-dropout LF. In or- 
der to investigate the evolution of the LBG LF over cos- 
mic time, we use the analyses of lSteidel et al.l (|1999l) and 
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iBouwens et alj (|2006D at redshift z = 3 and z = 6, re- 
spectively, since these are widely used for comparisons in 
other works, and represent an ideal compromise between 
area coverage and depth. These are plotted in Figure 
M The z 3 LF is K-corrected to 1400 A and adapted 
to our cosmology {M^,{z = 3) = —20.9); it is based on 
identical assumptions as in our fiducial LF, with which 
it can thus be directly compared. The z ~ 6 LF is based 
on a "cloning" algorithm, which is well bracketed by our 
two different assumptions for the dust extinction and the 
hydrogen intergalactic absorption. 

The possible underdensity of the HUDF, and the asso- 
ciated uncertainty in the faint-end slope of the z ~ 5 LF, 
is relevant in establishing the evolution of such faint-end 
through the comparison with the z ~ 3 — 6 estimates. In 
particular, the upper panel of Figure [H] shows that the 
inclusion of the HUDF in the computation of the z ^ 5 
LBGs LF leads to a shallower faint-end slope than at 
z ~ 3; in contrast, the faint-end slope that we derive us- 
ing the NICP12 field alone for the z ^ 5 LBGs is very 
similar to the one measured at that later epoch. A simi- 
lar result is found when the SMC extinction law and the 
Meiksin intergalactic extinction relation are assumed, as 
shown in the lower panel of the figure. This result is thus 
robust also towards possible variations of the dust ex- 
tinction properties through cosmic time. The possibility 
remains that the HUDF is, as we discuss above, ~ 30% 
under-dense relative to the NICP12 (and the mean); a 
renormalization of the HUDF by such an amount would 
actually lead to a slightly steeper faint-end slope at z ~ 5 
than at z 3. 

Rather similar conclusions are drawn for the redshift 
evolution of the faint-end of the LF in the z ~ 5 to 
z ~ 6 redshift window. Inspecting both panels of Figure 
[SI our fiducial z ~ 5 faint-end slope is consistently shal- 
lower than the z ~ 6 estimate, for both the HUDF and 
the NICP12 field; the NICP12 field provides however, at 
z ~ 5, a similar faint-end slope to the one measured at 
z ~ 6 when the SMC dust extinction law and the Meiksin 
intergalactic absorption correction are adopted. A sim- 
ilar agreement between the z ~ 5 and z ~ 6 faint-end 
slopes would also be obtained if the HUDF were renor- 
malized to account for a ~ 3 0% underdensity cause d by 
cosmic variance. Note that IBouwens et"all (|2006D did 
renormalize the z ~ 6 HUDF number density, increas- 
ing it by about 32%, based on the comparison with the 
GOODS data. 

8.2. Comparison with Theoretical Predictions 

We compare our observed z ~ 5 LBG LF with theo- 
retical predictions from the GALFORM semi-analytical 
model (jCole et al.ll2000f ). which calculates galaxy forma- 
tion in the framework of the CDM model of structure 
formation. It includes the assembly of dark matter halos 
by mergers, heating and cooling of gas in halos, star for- 
mation from cold gas, feedback from supernovae and pho- 
toionization, galaxy mergers, and chemical evolution of 
stars and gas. Star formation occurs in two modes: quies- 
cent in disks, and starbursts triggered by galaxy mergers. 
Galaxy luminosities are calculated based on stellar popu- 
lation synthesis, and include obscuration by dust, which 
is calculated self-consis tently by radiatiye transfer us- 
ing th e GRASIL model (|Silva et al.lll998l: iGranato et all 
I2OOOI) . The version of the model which we use is that 
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Fig. 9.— Evolution of the rest-frame 1400 A LBG LF in the 
redshift range 2 ~ 3 — 6. Top panel: The LFs assume the starburst 
dust extinction relation and the intergalactic absorption prescrip- 
tion of Madau (1995). Our fiducial 2 ~ 5 LF is represented by 
the black circles {NICP12 data) and the black squares (HUDF 
data), plus the black crosses, which represent the 2^5 SDF 
LF of Yoshida et al. (2006). In blue we show the 2 ~ 6 LF of 
IIBouwens et al.|[2006l) . and in red the 2 ~ 3 LF of llSteidel et al.l 
1 19991 V Solid and dotted lines are used to indicate the magnitude 
baseline over which the LF were measured and extrapolated, re- 
spectively. Bottom panel: A similar plot as above, in the assump- 
tion of a SMC dust extinction law and a Meiksin (2006) intergalac- 
tic absorption correction. The parameters for the LFs are listed in 
Table [3] 

described in iBaugh et al.1 (|2005l ) (see also iLacev et al.l 
l2007al ). and is based on the A-CDM cosmology, with 
f^o = 0.3, h = Ho/{100kms-^Mpc) = 0.7, and erg = 0.9. 
An important di fference of this model from the earlier 
I Cole et al] ()2000l ) model is the much large r role of star 
forma tion bursts at high redshift. In the iBaugh et al.1 
(|2005l ) model, bursts are responsible for 30% of the total 
star formation when integrated over all red shifts, com- 
pared with around 5% in the Col e et ahl model, and 
bursts dominate the total star formation rate density at 
z > 3. This difference in behaviour results from (i) mod- 
ifying the quiescent star formation time-scale in galactic 
discs, in order to make mergers at high redshift more 
gas-rich; and (ii) allowing triggering of bursts by both 
minor and major mergers. The predictions of this model 
fo r LBGs at d i fferent redshifts will be presented in full 
in ILacev et al.l ()2007bl ). here we give only the results for 
z ~ 5. The predictions of this same model for Lya- 
emitting galaxies at similar red shifts have already been 
given bv lLe Delliou et all ([2006). 

The Initial Mass Function (IMF) has a large impact on 
the predicted z 5 LF. We therefore computed the re- 
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suits for two different assumptions, n amely (1) a model 
where all stars are produced with a iKennicuttI (|1983l ) 
IMF, and (2) a model where quiescent star formation 
is described by a Kennicutt IMF, but star formation 
in the merger-induced starbursts occurs with a "Top- 
Heavy" (x = 0) IMF. The Kennicutt IMF has x = 0.4 
for m < Mq and x = 1.5 above this thres hold; it is 
therefore quite similar to a ISalpeterl ()1955[ ) IMF for 
m > Mq. An IMF similar to the Kennicutt IMF has 
been argued to better fit the observational data from our 
Gala xy and othe r spiral galaxies than a Salpeter IMF 
(e.g. IScalol 119981 ). At constant stellar mass, the far-UV 
(lOOO-SOOOA) luminosity for a Kennicutt IMF is within 
10-20% of that for a Salpet er IMF. The primary motiva- 
tion in iBaugh et al.l (|2005( ) for introducing a top-heavy 
IMF in starbursts was in order to explain the number 
of sub-mm galaxies seen at high redshift - with a nor- 
mal IMF, the predicted number counts of faint sub-mm 
galaxies were too low by more than an order of mag- 
nitude. However, a top-heavy IMF in st arbursts is als o 
motivated by theoretical calculations (e.g. lLarson II1998D . 
which suggest that, under the conditions that exist in 
starbursts, the IMF could be biased towards more mas- 
sive stars relative to quiescent star formation. 

Figure [10] shows the rest-frame 1500A L F at z ^ 5 
predic ted by the semi-analytical models (jLacev et al.l 
l2007bf l. compared to our measurements. The error bars 
on the theoretical curves show bootstrap estimates of the 
uncertainty due to the finite number of galaxies in the 
simulated samples; the numerical runs are very CPU in- 
tensive, and it is thus currently impractical to greatly re- 
duce these error bars by increasing the sample size. The 
dotted black and blue lines in the figure show, for the 
simple-Kennicutt and Kennicutt-|- Top-Heavy IMFs, re- 
spectively, the predicted LFs without dust obscuration. 
These substantially overpredict the number density of 
z ~ 5 LBGs over the whole luminosity range sampled 
by the data, but this is not surprising, since the average 
dust extinctions of star-forming galaxies in the far-UV 
at high redshift are thought to be quite large. 




Fig. 10. — Comparison of our fiducial z ^ 5 LF (red solid line: 
fit to the NICP12, HUDF and SDF data; red circles: NICP12 
data; red diamonds: HUDF data; black crosses: SDF data) with 
the prediction of Lacey et al. (2007). The simple-Kennicutt and 
Kennicutt-|- Top-Heavy IMF models are shown in black and blue, 
respectively; dotted and solid lines indicate unobscured and dust- 
extincted models. K-corrections from 1500 A to 1400 A are negligi- 
ble and were not adopted here. 



The solid black and blue lines are the simple-Kennicutt 
and Kennicutt-|- Top-Heavy IMFs models including the 
effects of dust extinction. At the bright-end of the z ~ 5 
LF, the simple-Kennicutt model with dust extinction 
substantially underestimate the number density of bright 
{Muv < -19) z - 5 LBGs. The Kennicutt-|-Top-Heavy 
IMF model is much closer to the observational data. Al- 
though this latter model still appears somewhat low com- 
pared to the observations, this could be explained by rel- 
atively small changes in the dust extinction, to which 
the model predictions are quite sensitive. Resolving this 
issue may require better observational estimates of the 
dust extinctions in LBGs. Another possibility would be 
that either the fraction of mass produced in bursts or the 
fraction of high mass stars in the bursts is even higher 
than in our Top-Heavy IMF model. 

The Kennicutt-|- Top-Heavy IMF model also appears 
to give a somewhat better fit to the faint end of the 
observed z ^ 5 LF than the simple-Kennicutt model, 
although the difference between the two models at the 
faint end is only slightly larger than the scatter in our 
estimate due to completeness corrections, input physical 
parameters, and possibly cosmic variance. 

9. SUMMARY AND CONCLUDING REMARKS 

We have presented the UDF05 project, a follow-up of 
the HUDF aimed at increasing the area of the HUDF 
and securing two more fields with ultra deep ACS and 
NICMOS imaging. 

In this first paper we use the UDF05 to estimate the 
z 5 Luminosity Function of Lyman-Break Galaxies 
down to very faint absolute magnitudes. Specifically, 
we use the NICP12 ACS field of the UDF05, and the 
ACS HUDF for a comparison, to extend the measure- 
ment of the LF of star forming galaxies at z ~ 5 down to 
Mi4oo = —17.1. This allows us to constrain the faint-end 
slope of the z ~ 5 LF very accurately to a ~ —1.6, which 
is in good agreement with the results of Giavalisco et al. 
(private communication) based on the GOODS data to 
a similar depth as in our work. After this paper was sub- 
mitted, Bouwcns et al. (2007) published LBG LFs from 
z ~ 4 — 6 based on all the available deep ACS fields and 
also their results are in very good agreement with our 
z - 5 LF. 

Attention was paid to correct the raw number counts 
for selection biases and photometric errors, and to esti- 
mate the impact on the LF of different assumptions for 
the properties of the z ~ 5 LBG population. We find that 
different assumptions about the LBG galaxy population 
and intergalactic absorption result in an uncertainty in 
the density of LBGs at z ~ 5 of about 25%. 

Under similar physical assumptions, the HUDF is un- 
derdense in z ~ 5 LBGs by about 30% with respect to 
the NICP12 field, a variation which can be accounted 
for by the expected amount of cosmic variance in the 
pencil-beam volumes we have probed. 

A substantial steepening of the faint-end slope is 
observed from th e local Universe out to z = 3 
(|Rvan et al.l [20071 ). in agreement with semi-analytical 
prediction s within a hierarchical galaxy formation sce- 
nario (K hochfar et a l. 2007). Further steepening of the 
faint-end slope a bove z '--^ 3 has been reported b y 
llwata et al.l (|2007l) . As discussed bvlOuchi et all (I2004D . 
however, the selection window of llwata et al.l ( 2OO70 is 
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rather broad, and it is thus hkely that their sample con- 
tains a lar ge n umber of low-redshift interlopers (see also 
Appendix IA.1[) . This would lead to an overestimation of 
the bright-end of the LF and thus to a spurious evolu- 
tion of the LF. Our much deep er data, exte nding about 
3 magnitudes fainter than the llwata et alJ (j2007i) sam- 
ple, show in contrast no evolution of the LF faint-end 
of LBGs in the z ^ 3 to 6 redshift window, and sub- 
stantiate a picture in which, once cosmic variance is 
taken into account, the faint-end of the LF of LBGs re- 
mains constant throughout this cosmic period. This is 
particularly interesting in comparison with the signifi- 
cant evolution shown, over the same redshift range, by 
the bright Lyman-break galaxy population (see also e.g. 



lYoshida et al1l2006HBouwens et al.ll2007f) . 

A comparison with our semi-analytical models, specif- 
ically aimed at investigating the effect of the IMF, sug- 
gests the possibility that the IMF of z ^ 5 LBGs is more 
top-heavy than a Salpeter IMF. Other factors could, 
however, be relevant, such as variations in the dust ex- 
tinction corrections. 
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APPENDIX 

A. COMPARISON WITH PUBLISHED Z - 5 LBGS LUMINOSITY FUNCTIONS 

A.l. The Subaru Deep Field 

Ground-based Subaru observations have been used to derive the z ~ 5 LF of LBGs (llwata et al ll2003l: lOuchi et al.l 
12004 lYoshida et al.l[2006l ). In particular. lYoshida et afl ()2006f ) extended the analysis of lOuchi et al.l ()2004h on the SDF 
with observations that reached 0.5 mag fainter and were wider in area by a factor of 1.5. As both authors adopted very 
similar techniques and their number counts of z ~ 5 galaxies are in very good agreement, we focus on the results by 
lYos hida ct al. (2006) only. Using the combination of Subaru V, i?, i' , z' filters, these authors defined two samples of 
z ~ 5 LBGs, whose reliability was tested with spectroscopic data. We can directly compare with their Ri'z' analysis, 
as this filter combination is sensitive to a very similar redshift wi ndow (z = 4.5 — 5.4) a s our data. 

Following a procedure similar to the one adopted in this paper, lYoshida et al.l ()2006l ) estimated the selection prob- 
ability and completeness by adding artificial galaxies to the original images; they did not correct their samples for 
photometric errors as we do, but estimated the number of expected low-z interlopers by s imulating galaxies with colors 
according to a photo-z catalog fr om the Hubble Deep Field North (jWilliams et al.iri996t ). 

The best- fitting LF derived bv lYoshida et al.l (|2006( l is listed in Table [3|(only Vi'z' sample, for the cases of keeping 
all LF parameters free and of fixing the faint-end slope to their z ~ 4 result). In both cases the faint-end slopes are 
very steep, i.e., significantly steeper than the one that we measure from our 3-magnitudes deeper data. This highlights 
the need for very deep imaging surveys to study the contribution to galaxy assembly at high redshifts by galactic 
systems well below the few bright perc ent of the star form ing population at those epochs. 

Note that we did not use the LF of llwata et al.l (|2003|) which has a larger amplitude due to a possibly significant 
contamination by low redshift interlopers (see also lOuchi et al.ll2004l : ICapak et al1l2004l ). 



A. 2. A Previous Analysis of the Hubble Ultra Deep Field 

iBeckwith et al.l (|2006f) present the first attempt to derive the -B435- to i775-dropout LFs in the HUDF. Their approach 
is very different from ours, as they did not correct for observational biases on the selection of dropout sources. 
Instead, they assumed that these biases affect all dropout samples equally, and base on this assumption their study 
of the evolution of the LBGs LF in the z ^ 4 — 7 redshift window. They extract an iyvs-selected sample from the 
GOODS vl data to constrain the bright-end of the LF, and derive the following best fit parameters: 0* = 0.001 
Mpc-^mag-i, M*(1400 A) = -20.5, a = -1.6 (fixed). 

Figure [TT] shows our fiducial z ~ 5 LF in comparison with the estimate of Beckwith et al. The two LFs are in good 
agreement, except for one datapoint at M = —18.2. This is a very interesting result, as it su ggests that the s evera l 
effects, which we have corrected for in our analysis, likely cancel out, thereby validating the IBeckwith et ahl (j2006l ) 
assumption, at least for the z ~ 5 LBG population. In a future paper (Oesch et al. 2007, in preparation) we discuss 
whether a similar conclusion can be drawn at z ~ 6. 



B. SIMULATIONS: DETAILS 

B.l. Idealized- Galaxies Simulations 

Two different types of light profiles were considered for the ideal-galaxies simulations, namely an exponential and 
a deVaucouleurs profile. These distributions are expected to more or less bracket the whole range of galactic light 
profiles. Specifically, the surface brightness distribution of disk galaxies was approximated by: 



ab 

where a and b are the major and minor axes lengths, the inclination angle i is given by cos(i) = ^, and Ltot is the 
extrapolated total galaxy light. The coordinates x and y are measured from the centre of the galaxy. The disk was 
assumed to be infinitely thin, which would lead to unrealistic shapes for edge-on systems. To avoid this unphysical 
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Fig. 11.— Comparison of our LF to the published HUDF LF bv lBeckwith et aLI I I2006I) . Our best fit to the combination of the HUDF 
data (red diamonds), NICP12 data (black squares), and the SDF data (blue crosses) is shown as a black line, and the estimates of Beckwith 
et al. are shown in green (symbols and line for the measurements and best fit Schechter function, respectively). 



situation, galaxy shapes were simulated with inclinations up to a maximum of 80°. To retain a flat distribution of 
inclinations in the entire 0° to 90° range, inclinations > 80°, drawn from an assumed flat distribution, were set equal 
to 80°. 

The deVaucouleurs profile adopted to simulate spheroidal-type galaxies reads: 



7.6692^ 



with r the radial distance as measured from the centre of the galaxy, and Ltot the extrapolated total galaxy light. 
Spherical symmetry was assumed for the deVaucoleurs galaxies. The r^/*-profile is very steep at small radii, and thus 
these types of galaxies are detected down to a fainter limiting magnitude with respect to disk galaxies. The total 
light is however distributed over a more extended region- and therefore, the amount of flux lost in the low-surface 
brightness wings is larger - than for an exponential profile. 

The idealized galaxies were convolved with the ACS PSF before being added to the NICP12 and HUDF images. The 
PSF was derived from a relatively bright but unsaturated star in the field. Tests with a PSF produced with TinyTim 
(|Krist. J.|[2004) . downgraded to the measured FWHM of the observed PSF, were also performed, and produced similar 
results. Poisson noise from the galaxy counts was not added to the simulated galaxies, since at the faint magnitudes 
of this study the sky background is the dominant source of Poisson noise. Tests were nonetheless performed, which 
confirmed that the galaxy shot noise had no impact on the results. 



B.2. Dimmed- Galaxies Simulations 

The idealized-galaxies simulations allow us to relate the intrinsic properties of galaxies to how these would be 
observed with our instrumental setup. They suffer however from the limitation of assuming perfectly smooth light 
distributions, as described by the adopted analytical descriptions. This is clearly not the case for real galaxies at any 
epoch, and certainly not for the UV light emerging from the star forming z ~ 5 galaxy population. 

In order to have a benchmark to gauge the incompleteness in our samples on the basis of more realistic assumptions 
concerning the galaxy light distributions, we performed an additional set of simulations, entirely based on the galaxies 
that were detected in the NICP12 and HUDF fields. Specifically, we dimmed the ACS images by specific amounts, 
and added these dimmed images to the corresponding original frames. Offsets of a few pixels were introduced before 
CO- adding the dimmed and original images. These offsets were chosen large enough to ensure that a given original 
galaxy was not 'self-blended' to its dimmed version in the subsequent SExtractor run. Since this procedure adds 
additional noise to the images, the co-added images were appropriately rescaled in order that galaxies of any given 
magnitude will have the same S/N values in the original and the final image. More specifically: 



1. A copy of the original image was divided by a factor /; 

2. This /-scaled version of the data was shifted by a few pixels relative to the original parent frame and added to 
the latter; 



3. The resulting composite image was divided by ^ 1 + (1//) , to retain the S/N vs. magnitude relation of the 
original image; 

4. The rms map of this final "dimmed image" was obtained from the rms map of the original image {(Jorig) and the 
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shifted image {u shift) by computing: 



^Ihiftl + '^orig 



The procedure leads to an effective dimming factor for the galaxies' fluxes of Udim — / ' y 1 + (1//)^ — ■\/l + 

We performed four separate experiments with four different shifts and dimming factors Udim — a/S, 2,4,6. For each 
experiment we shifted and dimmed the zgso, ^775, and Vgoe images by the same amount and reran SExtractor with the 
same parameters as for the compilation of the original catalogues. We used these " dimmed-galaxies" simulations to 
quantify, in an independent and complementary way relative to the " idealized-galaxies" simulations discussed above, 
the completeness of our catalogs, as well as the scatter in measured magnitudes and colors, as a function of galaxy 
luminosity. 

Note that these simulations differ conceptually from the 'cloning technique' adopted bv lBouwens "eFaLl (pOOeh as we 
add dimmed versions of the whole ACS images to the original tiles rather than inserting modified cut-outs of a subset 
of galaxies. 
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